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Abstract 

We explicitly compute the causal structure of the Schwarzschild black hole spacetime, by providing an 
algorithm to decide if any pair of events is causally related. The primary motivation for this study comes 
from discrete quantum gravity, in particular the causal set approach, in which the fundamental variables can 
be thought of as the causal ordering of randomly selected events in spacetime. This work opens the way to 
simulating non-conformally flat spacetimes within the causal set approach, which may allow one to study 
important questions such as black hole entropy and Hawking radiation on a full four dimensional causal set 
black hole. 
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1 Introduction 

Much is understood about the causal structure of the Schwarzschild black hole spacetime, e.g. in the sense 
that the collection of all null geodesies has been completely classified [IJ. However, what is missing is a com- 
plete specification of the causal relations, namely the set of ordered pairs of events in Schwarzschild which are 
connected by a causal curve. It is the purpose of this paper to provide exactly such a prescription: Given an 
arbitrary ordered pair of events in Schwarzschild spacetime, does there exist a future directed causal curve from the first 
to the second? 

Unfortunately, the differential equations describing the null geodesies are solvable in closed form only for 
limited values of their parameters HI. (In fact only a set of measure zero are expressible in closed form.) We 
therefore describe an algorithm which will allow one to compute an answer to the above question, to any (fi- 
nite) precision. For more implicit discussion on how the space of null geodesies encodes the causal structure of 
Schwarzschild, and general spacetimes, please refer to [2,^ and references therein. 

The primary motivation for this work comes from discrete quantum gravity, in particular the causal set ap- 
proach, for which the fundamental variables can be regarded as the causal ordering of randomly selected events 
in spacetime SUSllBj. More specifically the causal set is a discrete set of 'atoms of spacetime', which possesses 
a partial order relation which corresponds to causal ordering in spacetime. Because of this straightforward 
interpretation of the fundamental variables, it is relatively easy to extract phenomenological predictions from 
the theory on the effects of fundamental spacetime discreteness. Perhaps the most famous of these is the pre- 
diction of a small but non-zero fluctuating cosmological constant, whose current order of magnitude matches 
its observed value HZl- 

Mathematically a causal set is a set C endowed with an order relation -< with is irreflexive (x 7^ x), transitive 
(x -< y and y ~< z =^> x -< z), and locally finite (|{y|x -< y ~< z}\ is finite for all x, z G Cj^ The connection 
between a microscopic discrete causal set and a macroscopic approximating continuum spacetime arises via 
the notion of a 'sprinkling', which is a simple algorithm for generating a causal set from a spacetime. Given 
a spacetime M with finite spacetime volume V (such as a bounded region of an infinite spacetime), select at 
random events in M, with respect to the volume measure. By this we mean that in any region of M of 
volume V one expects to find n sprinkled events, where n is sampled from a Poisson distribution of mean v. 
Thus the probability of finding n sprinkled events in a region of volume v is 



(Therefore N is sampled from the distribution Eq. ([T]) with v = V .) Each of these sprinkled events then corre- 
sponds to an element of the causal set. After sprinkling the N events, one defines the partial ordering by stating 
that two elements are related iff the events are causally related in the spacetime. The microscopic-continuum 
correspondence then arises by the statement that if a causal set is likely to have arisen from a sprinkling into 
a given spacetime, then one regards that spacetime as a good macroscopic approximation of the underlying 
causal set. 

Thus far only conformally flat spacetimes, viz. those whose metric is given by a scalar conformal times the 
Minkowski metric, have been used in detailed calculations involving sprinklings of causal sets, because the 
second step of deducing the causal relations is particularly easy in that case. Here we open the possibility for 
sprinkling into a non-conformally flat spacetime: the Schwarzschild black hole. This allows one to address the 
general question of whether the methods of deducing properties of continuum spacetime from the causal set 

^The vertical brackets | ■ | stand for cardinality. 
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carry over to curved spacetime, and in particular spacetime with non-flat conformal structure. Some constructs 
which one may like to test are dimension estimators [5.8], the recovery of lengths of timelike 16^191 and spacelike 
UTOl geodesies, and extraction of macroscopic spatial topology |1JJ. 

The black hole is important for a number of other reasons, beyond merely being an example of a non- 
conformally flat spacetime. In the same way in which one can get a first approximation of the entropy of a 
gas merely by counting molecules, there is now substantial evidence that one can do an analogous counting 
of 'horizon molecules', such as causal links of a causal set crossing the horizon, to compute the entanglement 
entropy of a black hole, and thus get a handle on its physical origin Ill2| . In addition one may test a recently 
proposed entropy bound ||13| in this spherically symmetric, yet conformally curved context. Finally, this work 
can shed light on the longstanding problem of the role trans-Planckian modes in Hawking radiation, by allow- 
ing the study of a full four dimensional causal set black hole. 

The technique we present here, of deducing the causal relation from a classification of the null geodesies 
of a spacetime, can be generalized to other spacetimes as well. The most obvious example is the Reissner- 
Nordstrom black hole. It remains spherically symmetric, which considerably simplifies the analysis, and its 
null geodesies are completely classified [IJ. Presumably one could also analyze the Kerr black hole in this way, 
though the broken spherical symmetry will be considerably more complicated. Note also that one must avoid 
any regions which contain closed timelike curves (or consider a more general ordered substructure), since the 
causal set description breaks down there. 

In addition to describing an algorithm to decide if any two events in Schwarzschild are connected by a 
causal curve, we implement this prescription as a 'thorn' (module) in the Cactus high performance computing 
framework [14J. An advantage of doing so is that others can easily make use of the code, without having to 
write their own implementation of the algorithm, nor having to understand the details of the representation of 
the causal set on the computer. 

This paper is organized as follows. In section |2] we briefly describe the equations governing the null 
geodesies we will employ, and describe how to use them to determine if two events are causally related or 
not. In section |3] we present some of the causal relations as computed by our algorithm, and pictures which 
illustrate the causal sets which arise by sprinkling into a region of Schwarzschild spacetime. Section H] contains 
some concluding remarks. In Appendix |A] we give a proof that the null geodesies we integrate to determine 
the causal relations are those that arrive earliest. In Appendix |B] we describe the details of our implementation 
of this method, including how to sprinkle into the Schwarzschild geometry with uniform density. 

2 Null geodesies and the causal structure of a Schwarzschild black hole 

2.1 Preliminaries 

In this paper, we want to find a general recipe to determine unambiguously whether two events in four di- 
mensional Schwarzschild spacetime are causally related to each other. This problem has a simple answer in 
Minkowski spacetime, since it is straightforward to show that if — (^2 — * i ) ^ + (a^i — 2:2 ) ^ + (yi — J/2) ^ + (-^i — ^2) ^ < 0, 
two events Ei = {ti,xi,yi, zi) and E2 = {t2, X2, 1/2,22) are causally related, otherwise there is no causal relation 
between them. Things become more complicated in curved spacetime, where in principle we need to integrate 
the infinitesimal invariant distance ds along every possible path from one event to the other to see if there is a 
null or timelike curve (causal curve) connecting them. 

Nevertheless, it is possible to solve this problem, given an understanding of the collection of null geodesies 
in the spacetime. In the Schwarzchild metric, it is always possible to choose a time parameter that only increases 
towards future, i.e. (9f is everywhere a future pointing timelike vectorH Then given two events Ei and E2 with 

^ We will choose Eddington-Finkelstein coordinates, which have a such a time parameter, below. 
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Figure 1: Illustration of the timelike worldline 7, on which all spatial coordinates (including r) are constant, 
and two null geodesies from Ei to 7. 

ti < t2, the only possible causal relation between them is that Ei -< E2, which means that there is a future- 
directed causal curve from Ei to E2. Imagine a bunch of light rays (null geodesies) B, emanating from Ei, and 
E2 represents a particular moment t2 in the world line 7 of a stationary observer. If the world line of any null 
geodesies in B, which is a null geodesic, meets 7 at f < t2, then we can conclude that they are causally-related. 
On the other hand, if any null geodesic emanating from Ei can only reach 7 at t > t2, then Ei and E2 must be 
causally unrelated to each other. See Fig. [l] for an illustration. That this is true can be seen by proposition 2.20 
of [15J, which states if ^ G J {B) but A I^{B), then there must exist a future-directed null geodesic from 
Aio B. This tells us that the earliest future-directed causal curve must be a null geodesic, but since any null 
geodesic has failed to reach 7 early enough for a causal relation between the two events, then the conclusion is 
that there is no causal relation between them. 

We will prove in Appendix|A]that our procedure always considers the fastesjl geodesic from Ei to 7, so the 
arrival time t at 7 of that geodesic will be sufficient to determine if Ei and E2 are causally related. If we find 
that even this null geodesic meets 7 later than t2, we can say for sure that there is no way for these two events 
to be connected by any future directed causal curve. 

After this brief introduction, we discuss in the next subsection the particular simple case where Ei and E2 
are only radially separated, with no angular separation, so all we need to consider are radial geodesies. For 
more generic pairs of events, we must consider the full three dimensional casej^ which is too complicated for a 
complete analytic solution. However, before incorporating a numerical treatment for the generic case, in sub- 
section |23] we find two simple sufficient conditions for two events to be causally unrelated to each other. One 
uses a bound given by radial null geodesies, and the other by purely angular components. Furthermore, there 
is also a sufficient condition for two events to have causal relation, which is the existence of a composed null 

^By 'fastest' we mean the geodesic with the earliest arrival time, as given by the (EF) time coordinate. 

*Given two events in Schwarzschild spacetime, is is always possible to rotate the coordinates so that they both lie in the equatorial 
plane, as explained below. 
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curve connecting Ei and an event in 7 no later than E2. These sufficient conditions are enough to determine 
the causal relations for a large portion of pairs of events in Schwarzschild spacetime, and provide a very effi- 
cient preconditioning, since we only need to do numerical calculations for those which fail all these conditions. 
Furthermore, we shall give the recipe for generic pairs of events in subsection I2.4[ with the help of numerical 
calculations. As mentioned above, the proof of the lemma which ensures that the null geodesic found by our 
recipe is the fastest one connecting Ei and 7 is given in the Appendix. 

We start with the Schwarzschild metric in the familiar form 

ds^ = -(1 - ^)dtl + (1 - '^r'dr^ + r\dd^ + sin2 9dct?) , 

where M is the mass of the black hole. The metric is well known to possess a coordinate singularity at r = 2M, 
where the event horizon lies. It can be written in Eddington-Finkelstein (EF) coordinates [16], 

ds'^ = -dt^ + dr^ + r^(d9^ + sin^ edcp"^) + —(dt + drf , (2) 

r 

with the following transformation of time parameter. 



r 



t = ts + 2Mlii( 1) . (3) 

Given two events in the EF coordinates, Ei = (ti, ri, 0i, and E2 = (^2, ^^2, ^'2, ^2) with ti < t2, the 
only possible causal relation is Ei -< E2. Besides, it is obvious that one can always choose suitable angular 
coordinates and ip for which ■!?! = '!?2 = 7r/2, ipi = 0, and ip2 = arccos(cos 9i cos 62 + sin Oi sin O2 cos((/>i — ^2)) G 
[0, vr]. Therefore, it is sufficient to consider a pair of events Ei = (ti, ri, 7r/2, 0) and E2 = (^2, ^"2, vr/2, ip2 G [0, vr]) 
with ti < t2. The stationary worldline 7 containing E2 is given by (t,r = r2,'& = ^,(p = 'P2), i-e. all spatial 
coordinates are held fixed. We define the angle through which our null geodesic travels, /S.ip, by 



A(/7 = dip , 



V 



where is a null geodesic from Ei to 7, and ip is the azimuthal coordinate in the rotated coordinate system. 
In subsection l2.4l we shall find that a generic null geodesic from Ei to 7 can always be formulated as one with 
= 7r/2 as well, thus the whole problem, including the events and geodesies in between, can be projected to 
2 + 1 dimensions. 

An important question when considering the generic pair of events will be if the null geodesic we choose is 
the fastest one, i.e. that it arrives at 7 before all other null geodesies from Ei. In order for the null geodesic r] to 
arrive at 7, it must have Aip = Ikir + ip2, with A; G Z. In Appendix lAl we prove that a geodesic r] with \Aip\ = ip2 
is a fastest future directed null geodesic from Ei to 7!! Thus we only need to consider those null geodesies from Ei 
to 7 which travel for an angle no more than vr in the ip direction. 



2.2 Radially separated pairs and radial null geodesies 

Now let us consider the simplest case with ip2 = 0. In this case we only need to consider radial null geodesies, 
and it is straightforward that in the EF coordinates, by setting d-d = dip = and ds^ = 0, radial null geodesies 
take a simple form, 

dt + dr = 0, (4) 
/2M \ , /2M \ 

i — -ljdt+i — + ljdr = 0, (5) 

'The absolute value on \Aip\ is merely to account for the case of ip2 — 11,111 which there are two fastest null geodesies, i.e. that both 
have the same arrival time. However, note that this case will never arise in the sprinkling described in Section [l] as the set of such 
sprinklings is of measure zero in the space of all sprinklings. 
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where the first one describes ingoing null geodesies, while the second one turns out to be outgoing for null 
geodesies outside the event horizon, and ingoing for those inside the event horizon. The two radial directions 
for null geodesies, both outside and inside the horizon, are illustrated in Fig. |2l Therefore, given two events 



Singularity 



Horizon 





Exterior Region 



r 



Figure 2: Schematic illustration of lightcones in a Schwarzschild black hole. The event horizon is indicated by 
the blue line, and the interior region is to the left. Note the two ingoing radial directions for future directed 
lightcones in the interior region. (The vertical direction is to be regarded here as Schwarzschild time, which 
can go backward for timelike observers inside the horizon.) 

without angular separation, Ei = (ti, ri, 7r/2, 0) and E2 = (^2, ?'2, 7r/2, 0) with ti < t2, we can determine if they 
are causally related by considering the time for null geodesies from Ei to reach a point in 7, the world line of a 
stationary observer at (r2, 7r/2, 0). 

First assume ri > r2. We consider only ingoing null geodesies from ri to r2. If both events are outside the 
horizon, ri > r2 > 2M, there are only ingoing null geodesies given by Eq. and it is straightforward to see 
that these two events are causally related to each other if and only if t2 > ti + ri — r2, where 1 1 + ri — r2 is simply 
the time when a null geodesic from Ei hits 7. If ri > 2M > r2, we still only have Eq. @ for ingoing rays outside 
the horizon, and to ensure continuity of the tangent of the geodesic when it crosses the horizon, we can only use 
Eq. lH) inside the horizon too, and the conclusion is not altered. However, if both events are inside the horizon, 
2M > ri > r2, we have two choices of ingoing rays given by Eq. (H) and Eq. (|5]), and a simple calculation shows 
that it takes a shorter time for a null geodesic from Ei to hit 7 for the first case, dt = —dr = \dr\, than for the 
second case dt = — 2m^/1'^\ dr = (1 + 2M-r )\'^'''\ > \dr \ for r < 2M. Therefore, for both events inside the horizon, 
the necessary and sufficient condition for their causal relation is still t2 > ti + ri — r2. 

On the other hand, if r2 > ri, we shall instead consider outgoing rays, given by Eq. (O outside the horizon. 
There is a special ease, ri < 2M, for which no matter r2 < 2M or r2 > 2M, they can not be causally related to 
each other for any ti and t2, since there are no outgoing null geodesies either propagating inside the horizon or 
coming out of it. If they both lie outside the horizon, r2 > ri > 2M , then the sufficient and necessary condition 
for a causal relation between them is ^2 > + ^"2 - ri + AM ln(^j5|^), where h + r2 - n + AM ln( ^^I^*^ ) = 
ti + J^^ r-2M time for a null geodesic given by Eq. (O to travel from Ei to the point it meets 7. 

To summarize, we have the following recipe to determine if two events are causally related to each other. If 
i"i > r2, the (necessary and sufficient) condition for a causal relation is always 



t2>ti + ri-r2, 
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2.3 Sufficient conditions for causally related and unrelated pairs 

Since in practice we are interested in computing the causal relations between many pairs of events, it is useful to 
derive bounds that allow one to decide whether the events are related quickly based on simple criteria, without 
having to perform any numerical integrations. In this section we derive such bounds, which are sufficient to 
determine that the two events are unrelated (section [2.3.1|) or related (section |2.3.2| |. 



2,3.1 Spacelike bounds 

For two events with non-zero angular distance, i.e. 1^2 S (0, vr], we have to consider non-radial, three dimen- 
sional null geodesies in general [T|. Nevertheless, we shall find two independent sufficient conditions for the 
events to be causally unrelated below, by using lower bounds on the time duration of any future-directed 
causal curve from Ei to 7, given by either radial null geodesies with the angular part discarded, or by a purely 
angular component with the radial part discarded. 

Since anywhere on a future-directed causal curve we have ds^ < and dt > 0, with Eq. ^ we obtain an 
inequality by discarding the (positive) angular part. 



{dt + dr) 



2M \ f2M 
l\dt+[ + l\dr 



r I \ r 



< 



If r > 2M, we obtain for an ingoing causal curve dr < 0, dt > —dr, and for an outgoing curve dr > 0, 
dt > ^^2M C?r. If r < 2M, we obtain for an ingoing curve dr < 0, y^^dr > dt > —dr, and no solution for an 
outgoing curve dr > 0. 

These inequalities can be integrated along the whole curve from Ei to 7 to give a lower bound on the time 
duration, which turns out to be exactly the time duration of radial null geodesies from ri to r2. If ri > r2, the 
sufficient condition for two events to be causally unrelated is 

t2-ti <ri-r2 , (6) 

while if r2 > ri > 2M, they are unrelated if 

ro - 2M 

t2-ti<r2-ri+m\n{- — -) . (7) 

ri — IM 

Finally if r2 > ri and ri < 2M, they can not be causally related to each other. 

On the other hand, we can also try to discard the radial part to obtain a bound by using a purely angular 
component, but this is more subtle and we need to be careful. Since dr"^ + ^ {dt + dr)"^ > 0, we can discard this 
part to obtain an inequality anywhere on a causal curve, as —dt^ + r^d93^ < 0, which gives > r|(i(^|. However, 
for outgoing causal curves outside the horizon, we have (1 + ^)dr'^ + ^dtdr > and a stronger bound can 
be obtained as - l)dt'^ + r^dip"^ < which gives dt > -^=^=, where r > 2M. The equality holds for null 
geodesies with constant r, which satisfy 



1 - —dt ±rd(p = , (8) 
r 

where r = const > 2M. This can be rephrased physically by stating that it is impossible for a particle inside 
the horizon to move around an orbit with constant r < 2M, because nothing can stop it from falling into the 
singularity r = 0. We shall make direct use of Eq. ^ later for the timelike bound. 

For any ingoing curve from Ei to 7, we have dt > r\dip\ > r2\dip\ along it, and we obtain a sufficient 
condition for two events with ri > r2 to be causally unrelated, 

t2 <ti+ r2f2 ■ 
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For any outgoing curves from Ei to 7, we have dt > ^j^^'^a// ^^°^S need to find out the minimum 

of f(r) = , ''' in the range 2M < ri < r < ro. It is straightforward to obtain 

fir)- 



[1 - 2M/r)3/2 ' 

from which the location of the minimum can be determined to be ro = ri for 3M < ri < r2, tq = 3M for 
ri < 3M < r2, and rg = r2 for ri < r2 < 3M. 

Therefore we have dt > , ''"1'^'^^ along the curve and we obtain a sufficient condition for two events with 

- ^l-2A//ro ^ 

2M < ri < r2 to be causally unrelated, 

t2 - ti < /(ro)v'2 . (9) 

2.3.2 Timelike bound 

Furthermore, it turns out that for many pairs, as long as at least one of them is outside the horizon, we can use 
radial null geodesies and null geodesies with constant r to find a composed null curve connecting E\ and an 
event in 7. If this event is no later than E2, then this will be a sufficient condition for their causal relation. 

Given Eq. (H)), Eq. (O and Eq. (|8]), we can construct a null curve from E\ to an event in 7 whenever r\ > 2M, 
which is composed of a sequence of null geodesies with constant and constant r. To optimize this sufficient 
condition for a causal relation, we need to minimize the time duration of the segment of the null geodesic with 
constant r, which occurs at r = ro where ro = min(ri,r2) for ri,r2 > 3M, rg = 3M for r\ > 3M > r2 or 
r2 > 3M > r\, and ro = max(ri, r2) for ri, r2 < 3M. 

If ri > ro > r2, then we can compose the null curve by the following three segments: 

1. an ingoing radial segment from r = ri to r = ro with (/? = 0; 

2. a segment from 99 = to 99 = (/?2 with r = ro; 

3. an ingoing radial segment from r = ro to r = r2 with Lp = ip2- 

The time for this null curve to reach 7 is easy to compute by Eq. (|4| and Eq. (|8]), 

t = ii + n - r2 + /(ro)v32 • (10) 

Here we do not care about whether r2 is inside or outside the horizon, since regardless of whether the segment 
3 lies completely outside the horizon or crosses the horizon, we always use Eq. (H)) for continuity. Similarly, if 
^^2 > ''o ^ ^'1 > 2M, then we can construct the null curve by replacing segments 1 and 3 above by outgoing 
radial segments, and the time can be computed by Eq. (O and Eq. ((8]l to be 

r9 - 2M 

i = ti + r-2 - ri + 4M ln(^ — ) + /(ro)v92 • 

r\ — 2M 

Therefore, we have the following sufficient condition for two events, of which at least one is outside the 
horizon, to be causally related. If ri > r2 and ri > 2M, then they are causally related if 

t2>h+ri-r2 + f{ro)(p2 ; 
if r2 > ri > 2M, they are causally related if 

ro - 2M 

t2>ti + r2-n+ 4Mln{- — -) + f{ro)ip2 . 

ri — 2M 

If a pair of events fails both of these sufficient conditions, then we have to consider a generic form of null 
geodesies to determine if they are causally related. 
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2.4 Generic pairs of events and null geodesies 

The most generic null geodesies in Schwarzschild spacetime have the following form 

dt dr d6 dcj) 

Pt-; Pr-; Pe-j P,t>-r = ' 

dr dr or ar 

where r is an affine parameter, and pt, Pcj) are constants, 

2M dts p 

Pt = [l H- = E, 

r dr 

p^ = r sm 0— = L , 

and p^ satisfies 

-^^^ = r sm d cos 0[ — ) . 

dr dr 

If we choose -d = 7r/2 at the moment when ^ = 0, we get also ^ = at this moment, which implies ?? = 7r/2 
all along the geodesic. Therefore a general null geodesic can be described in the plane ^9 = 7r/2, which also 
simplifies its equations to 

(*,=+^(l-^, = E^ («) 

ar r 

where E and L denote the constant energy and angular momentum of the massless particle, 

r dr 



dip L 
dr r 



(12) 



The energy is expressed in terms of the Schwarzschild time parameter ts- 

The full set of Eq. ((TT]| and Eq. ((T2ll , combined with initial values of t, r and (p, as well as E and L, can 
uniquely determine a null geodesic in Schwarzschild spacetime. However, since we only want to obtain rela- 
tions between t, r and ip without the affine parameter r, it is convenient to consider r as a function of ip and use 
a new variable u = 1/r. We then obtain from Eq. (|TT]) , 



or equivalently. 



dip 



^ = ±{2Mu'-u'' + c')-'/\ (13) 
du 

where + corresponds to dip/du > 0, — corresponds to dip/du < 0, and c = E/L for L 7^ 0. The case of L = 0, 
which corresponds to radial null geodesies, has been discussed in Subsection 12.21 It turns out that, for L 7^ 0, 
the geodesic depends on E and L only through their ratio c. Using Eq. lO and Eq. ((12)), we can further obtain 

dt cr"^ ^ dr/dip 



dip l-2M/r r/2M-l ' 
which can be simplified by using the new variable u = l/r and Eq. ([131) as 



_ c =F 2MuV2Mu^ - + 
d^~ V? - 2Mu^ ■ ^ ' 
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Alternatively, we can put it into an equation involving only t and u, 



dt _ ±c{2Mu^ -u^ + c2)-i/2 _ 2Mu 

d^ ~ n2 - 2Mn3 ' ^ 

where, as mentioned before, + corresponds to dip/du > 0, — corresponds to dip/du < 0. Eq. ((T3)l and Eq. ((141) . 
or equivalently Eq. ([T3l > and Eq. dlSl ) is a full set of equations for generic null geodesies with non-zero angular 
momenta. 

Now given Ei and E2 which do not satisfy any of the sufficient conditions in Subsection |23J we have to do 
the following numerical calculation to see if they are causally related to each other. For any c we can integrate 
Eq. ((T3|) from ipi = 0, ui = 1/ri to U2 = l/r2, and get some value (^2. By choosing a suitable cq, we can make 
^'2 = ^2 which means that the null geodesic with cq hits 7 from Ei . Then we can use cq in Eq. ((T5)l , and integrate 
it from ui, ti to U2 and get some value If t < t2, then they are definitely causally related. If i > t2, according 
to the lemma of Section l2ll they must be causally imrelated. 



3 Results 

3.1 Causal Relations 

We sprinkle into a region of Schwarzschild spacetime, which is bounded by < r < rmax arid < t < tmax- See 
Appendix IB . 1 1 for details on how this is done. It is important to note that in this section we use the unrotated 
coordinates, so 9 is not restricted to |, and use coordinates such that M = 1. By 'equatorial plane', we simply 
mean that 6* = f • 

Nine events selected from a region of Schwarzschild with rmax = 3, tmax = 8, are shown in Table [10 



event 


t 


r 


e 


4> 







0.410895 


2.36161 


1.80295 


0.57951 


exterior 


1 


1.109415 


2.89891 


1.04335 


4.25531 


exterior 


2 


1.133105 


1.36083 


1.89919 


1.06482 


interior 


3 


2.743428 


2.74093 


2.97906 


4.22204 


exterior 


4 


3.235970 


0.65462 


0.11664 


5.06884 


interior 


5 


3.972871 


0.96354 


2.33727 


1.38169 


interior 


6 


5.230757 


2.34476 


1.11855 


3.47242 


exterior 


7 


6.014261 


0.664739 


2.82235 


0.95459 


interior 


8 


6.193089 


0.429636 


2.20122 


1.99644 


interior 



Table 1: Nine events in Schwarzschild spacetime, specified by their Eddington-Finkelstein coordinates. 

Our task is to decide, for each pair of events {Ei, E2) (with Ei having an earlier EE time coordinate than 
E2), whether they are causally related or not. To do this we perform the following algorithm: 

1. Is £'1 is behind the horizon and r2 > ri? If so they are unrelated. 

2. Change the angular coordinates so that both lie on the equatorial plane, and restrict attention to null 
geodesies which traverse an azimuthal angle < vr on their trip from Ei to a stationary worldline 7 con- 
taining £'2- 

3. Is the EE time separation of the events less than the angular or radial spacelike bounds? If so they are 
unrelated. 

*We chose rmax = 3 to get roughly the same number of interior and exterior events, and tmax = 8 to be roughly half the circumfer- 
ence of a circle at r = 3. Note that r = 3 corresponds to the innermost circular orbit (which is therefore lightlike). 
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4. Is the EF time separation greater than the timelike bound? If so they are related. 

5. If neither sufficient condition is satisfied, then we must numerically compute the value of {E/L)"^ which 
will send a null geodesic from Ei to 7. Armed with this value, we compute the elapsed coordinate time 
along this geodesic, and decide if it arrives before or after the event E2. 

Table |2] gives details for a selection of pairs of events from Table [T] In particular we show the various 
quantities which are computed along the way to deciding if this pair is causally related. For each pair, if 



pair 


dir 








V2 




At 


rad trip 


ang bnd 


tot trip 


{E/Lf 


time 


result 


01 


out 


2 


898906 


2 


567258 





698520 


4.179694 


13.36484 


— 


— 


— 


unrelated : either bound 


02 


in 




— 





475629 





722210 


1.000779 


0.647253 


— 


— 


— 


unrelated : radial bound 


12 


in 




— 


2 


937662 





023690 


1.538071 


3.997673 


— 


— 


— 


unrelated : either bound 


04 


in 


2 


361614 


1 


827161 


2 


825075 


1.706999 


1.196088 


12.73426 


0.0460462 


4.69799 


generic, crossing, unrelated 


14 


in 


2 


898906 





965528 


2 


126556 


2.244290 


0.632050 


— 


— 


— 


unrelated : radial bound 


24 


in 


1 


360835 


1 


973603 


2 


102865 


0.706219 


1.291951 








generic, hits singularity 


05 


in 


2 


361614 





867289 


3 


561976 


1.398076 


0.835666 


6.632326 


0.272754 


2.11915 


generic, crossing, related 


15 


in 


2 


898906 


2 


821962 


2 


863456 


1.935368 


2.719068 


16.62616 


0.0388018 


9.86353 


generic, crossing, unrelated 


25 


in 


1 


360835 





512295 


2 


839766 


0.397297 


0.493616 








generic, hits singularity 


45 


out 






















emerg 


ing from interior : unrelated 


06 


in 


2 


361614 


2 


820685 


4 


819862 


0.016859 


6.613817 








unrelated : angular bound 


16 


in 


2 


898906 





690536 


4 


121342 


0.554150 


1.619139 


4.148999 


0.0476468 


2.60973 


generic, exterior, related 


26 


out 






















emerg 


ing from interior : unrelated 


46 


out 






















emerg 


ing from interior : unrelated 


56 


out 






















emerg 


ing from interior : unrelated 


37 


in 


2 


740935 





480906 


3 


270833 


2.076196 


0.319677 


4.611431 


4.55373 


2.1833 


generic, crossing, related 


58 


in 





963538 





485959 


2 


220218 


0.533902 


0.208785 




1.35646 


0.667059 


generic, interior, related 



Table 2: Considering relations between a collection of pairs of elements of Table [H 



applicable, we show: 

• the direction (ingoing or outgoing) 

• the angle between the events 922 

• the time coordinate separation between the events At 

• the 'radial' spacelike bounds, given by Eq. ^ for ingoing null geodesies, and Eq. (O for outgoing geodesies 

• the 'angular' spacelike bound, Eq. ((H) 

• the coordinate time traversed along the trip, composed of successive segments at constant r or angular 
position, which yields the timelike bound of Subsection l2.3.2l 

• the value of (? = {EjVf' in Eq. ((13)) which yields the fastest null geodesic from E\ to 7. These values are 
computed using the procedure detailed in Appendix |B.2[ 

• the coordinate time elapsed along this fastest geodesic 

• whether they are related or not, and which condition allows us to decide 

For the 1 pair, the time coordinate separation At is less than either of the spacelike bounds, so the events are 
unrelated. The 2 pair fails the angular spacelike bound, but passes the radial bound. The 4 pair is 'generic', 
meaning that it fails both spacelike bounds and the timelike bounds. These simple tricks are insufficient to 
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determine if the events are related, so we must integrate Eq. (T3[ to locate the fastest null geodesic from event 
to the 7 containing event 4. This geodesic crosses the horizon, but does not arrive at 7 in time for the events to 
be related (since the 'time' in the last column is larger than the 'available time' At). The 2 4 pair fails the radial 
spacelike bound. The radial and timelike bounds fail because the pair is inside the horizon, for which there is 
no timelike, constant-r trajectory. Integrating ((T3|) . with = {E/L)"^ = 0, gives only (p = 0.721811, which is 
not enough to reach 7. All null geodesies will hit the singularity before reaching the 6*2 = 0.11664, 4)2 = 0.11664 
worldline. The 2 5 pair meets a similar fate: there are no future directed null geodesies from event 2 which 
reach 62 = 2.33727, 02 = 1.38169 before falling into the singularity. The 5 pair fails all bounds, and thus is 
generic. This time, however, the null geodesic does reach 7 before event 5. The 4 5 pair represents an attempt 
to 'escape from the interior', in the sense that event 4 is inside the horizon, and event 5 is at a larger radius than 
4. Even though event 5 is also inside the horizon, there are no causal curves inside the horizon which extend 
to larger radii. Events and 6 are at almost the same radius, but at very different angular positions. Thus the 
angular bound is useful in deducing that they are unrelated, without having to integrate any geodesies. Events 
1 and 6 are an example of a generic pair which are both outside of the horizon. They happen to be related. 
The pairs 2 6, 4 6, and 5 6 suffer the same fate as 4 5. This time they are even attempting to escape across the 
horizon. The 5 8 pair is generic and completely inside the horizon. This time there are causal curves which 
reach 7 from event 5, and the events end up being related. 



3.2 Causal Sets 

In this subsection we show some Hasse diagrams of the causal sets which arise from sprinkling into Schwarzschild. 
In the Hasse diagram one shows only the links of the causal set, namely those causal relations which are not 
implied by transitivity. Figures |3] and |4] use the graphviz package [17J to generate the diagram, ignoring the 
embedding information. Figures |5l HI [71 and [8] arrange the causal set elements using their embedded location. 




Figure 3: Hasse diagram of a 38 element causal set sprinkled into Schwarzschild spacetime, with rmax = 
4, ^max = 8. The circled elements landed behind the horizon. Note that no information escapes from this 
interior region, in that there are no causal relations from circled to uncircled elements. 



0000 ©00 (iZT © O® 0©0©® ©O©® ® ©® ® 

Figure 4: A 'causal set singularity'. This 64 element causet is sprinkled near the singularity of Schwarzschild, 
?"max = 1, imax ^ oo. (The resulting causet is independent of tmax beyond some finite value.) 

Figure |3] shows a 38 element causal set, which arises from a sprinkling into a region of a Schwarzschild 
spacetime. Figure H] portrays a causal set which arises when sprinkling near the singularity. Note that there 
is no need to worry about sprinkling on top of the singularity, as that is a zero probability event, even if the 
singularity is contained within the sprinkling region. It is very antichain-like, as the futures of the elements 
rapidly fall into the singularity, so it is unlikely that another sprinkled element lands in that region. 

In the remaining four figures of this section we show more Hasse diagrams, though this time we use the 
embedding information to locate the nodes in the graph. For the first three. Figures |5l |6l and [7] we generate a 
causal set C by sprinkling into half of the equatorial plane (0 = 7r/2, G [0, vr]), with tmax = 8 and rmax = 4. The 
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Figure 5: Hasse diagram of a causal set C which arises from sprinkling 91 elements into the half equatorial 
plane. Elements sprinkled behind the horizon appear blue, while those outside are green. 



13 




Figure 6: The same 91 element causal set C as viewed from the 'top' (t — > +00, look into the past direction). 
The blue 'half disk' corresponds to the black hole interior. 
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Figure 7: The sprinkled 91 element causal set C, this time plotting r against (j) interpreted as Cartesian co- 
ordinates, rather than polar coordinates. The blue region on the left corresponds to the interior of the black 
hole. 
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relations (links) are shown in red. The half plane is chosen to reduce the 'clutter' from links between distant 
elements (such as links which cross the black hole, from the exterior region on one side to the exterior region on 
the other). Note that these red lines are not null geodesies connecting the elements (though they will converge 
to such in the infinite sprinkling density limit). They are simply straight lines between the elements drawn by 
the plotting program (gnuplot). 

The last figure [8] attempts to more clearly show the light cone structure of Schwarzschild, in particular the 
tilting of the light cones toward the singularity at r = 0. Here we sprinkle 91 elements into a 1/10 radian 
'wedge' of the equatorial plane. Since the number of elements sprinkled into a region grows quadratically with 
r, we plot here t versus r^, so the distribution will appear approximately uniform in the horizontal direction. 
On the right side of the figure we have a cutoff at r = 4, so obviously no links can go beyond that cutoff. 
However around r = 3 we can see the light cone spreading in both directions. As we get closer to the center 
the links extend less and less to the right (exterior). As we cross the horizon at r = 2 we see that the links only 
extend to the left, since their are no future directed causal curves which have non-decreasing r in the interior. 




Figure 8: Sprinkling of 91 elements into a 1/10 radian wedge of the equatorial plane. Here we ignore the (/> 
coordinate, and square the radial coordinate in order to spread out the elements at large radius. 
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4 Conclusion and discussion 



We have described an algorithm to determine if two events are causally related in Schwarzschild spacetime. 
It involves first checking a number of sufficient conditions, to see if these are able to determine whether they 
are related, without resorting to a time consuming numerical integration. If none of these are satisfied then we 
locate a future directed null geodesic which leads from the earlier event to a worldline containing the later, and 
integrate the elapsed EF time along the geodesic to determine if the pair is related. 

We then wrote a 'thorn' (module) within the Cactus framework with sprinkles events into a region of 
Schwarzschild, and use this algorithm to deduce the causal relations between every pair of events. This proce- 
dure yields a causal set which 'faithfully embeds' into a Schwarzschild black hole spacetime. 

It should not be difficult to generalize this prescription to other black hole spacetimes, such as Reissner- 
Nordstrom or Kerr. In order to describe the closed timelike curves in interior regions one may generalize the 
definition of the causal set slightly, replacing the irreflexive condition with reflexivity, so one ends up with 
a transitive directed graph. In addition one can consider other conformally curved spacetimes, such as the 
/c = ±1 Friedman-Robertson- Walker universes. 

This work is a good starting point for one to address a wide array of questions within the causal set program, 
which previously were inaccessible. On one hand it will allow one to easily investigate many kinematical 
questions with regard to the so-called 'Hauptvermutung' of causal sets, which conjectures that if one has a 
causal set which is likely to arise from sprinkling into two separate spacetimes, then those spacetimes must 
be approximately isometric. To date, all results regarding to how to deduce properties of an approximating 
continuum from a causal set consider only conformally flat geometries. Now one will be able to test such 
constructs on a much wider class of geometries. 

Another important application is towards our understanding of black hole thermodynamics. The funda- 
mental discreteness of causal sets gives us a possible way to characterize the degrees of freedom which give 
rise to black hole entropy. We are now in a position to repeat the analysis of the link counting and its gen- 
eralizations llT2l , for the full 4d Schwarzschild geometry. As argued by Dou and Sorkin, this may provide 
access to the fundamental length scale of quantum gravity, since entropy, as a pure number, is not subject to 
renormalization. 

Finally, a long standing question in semi-classical gravity is the trans-Planckian problem, that the Hawking 
radiation emanating from a black hole horizon, at late times (after the black hole forms), arises from modes of 
frequency much greater than the Planck scale ITSl . If spacetime is discrete at this scale, then one may expect 
that such trans-Planckian modes cannot exist. How then is Hawking radiation possible, in such a discrete 
setting? Important groundwork on the dynamics of scalar fields on a background causal set has recently been 
laid l,19i . .20..21| . Now that we can construct causal sets which correspond to a full four dimensional black hole, 
it may be possible to address this question within the causal set approach. 
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A Appendix: Proof of the proposition 



In this section, we want to prove the following proposition. Given two events Ei and E2 in Schwarzschild 
spacetime, with ti < t2, and E2 representing the moment t2 on the world line 7 of a stationary observer, the 
null geodesic with \Aip\ = (p2 is the fastest one (arrives at the earliest time in EE coordinates) for all future- 
directed null geodesies from Ei to 7. 

It is clear that the fastest geodesic will have the least elapsed time. Thus from Eq. |(15)). we wish to minimize 
the integral 

along the geodesic. Since dtp is clearly nonnegative, the sign in front of c will be positive for ingoing geodesies 
{du > 0) and negative for outgoing geodesies (du < 0). Given Ei and E2 (and thus 7), different null geodesies 
are given by different values of c, and it is straightforward to see from Eq. ((T6)l that in every case for larger c. 
At gets smaller. Therefore, the fastest null geodesic from Ei to 7 must have the largest possible c to reach 7, 
which, by Eq. ([13]), has the smallest possible \ Aip\ = ip2^ 



B Numerical Details 

B.l Sprinkling into Schwarzschild Spacetime 

In implementing these ideas on a computer, the first task is to randomly select events in spacetime, with a 
density proportional to the spacetime volume factor 

= r-\/r2 + 12M2 sin 9 . (17) 

Because of the simple product form of this expression, we can break up the sprinkling into an angular piece, 
a temporal piece, and a radial piece. For the uniform sprinkling of the angular coordinates, on a 2-sphere, we 
follow the procedure described in Section 5.2 of 1.22,1 . For the temporal sprinkling, since the volume element is 
independent of t, we can select its values uniformly at random. 

The sprinkling in the radial direction must be performed such that it yields the distribution of Eq. ((17|| 
(ignoring the 9 dependence; this is accounted for above). This is achieved by the following general method 
(derived from [23J). To sprinkle a coordinate x between the bounds a < x < b such that it has distribution f{x), 
compute the indefinite integral 

I{x) = ^ iy{x')dx' , (18) 

where the constant normalization factor is 




Now invert Eq. flSl) to get x as a function of /. This expression, where I is a random variable distributed 
uniformly in the unit interval [0, 1], will be distributed according to f{x). Of course this method only works 
if f{x) is integrable and its integral is invertible. In our case /(r) = r-v/r^ + 12M2, and the expression is 

V(3iV/ + (a2 + 12Af 2)3/2)2/3 _ 12Af2.|8| 

^As mentioned earlier, for (^2 < ti", there is only one fastest null geodesic with the largest c to make A(p = (^2; while for <^2 ~ tt, 
there are in fact two fastest null geodesies with A(p — ±tv but the same c and At. 

^In fact the angular sprinkling is of this type as well. With f{6) — sin 6 the expression is arccos(2/ — 1). The (j) coordinate is 
distributed uniformly in [0, 27r]. 
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B.2 Determining the Causal Relations 

Once equipped with a collection of events in Schwarzschild, we can sort them by their time coordinate, and then 
consider each sorted pair in turn. The first task is to check the sufficient conditions described in Subsection |2.3[ 
This is relatively straightforward; numerous examples are given in Subsection 13. II If the pair fails all available 
sufficient conditions, then it is a 'generic pair', and we must integrate null geodesies as described in Subsection 



The basic task is to find a value of the parameter = {E/L)'^ for which the integral of Eq. ((13)) along the null 
geodesic from ui to U2 equals 992- This is made complicated by the fact that the cubic f{u) = 2Mu'^ — + in 
the denominator of the right hand side can have real roots within the domain of integration. (See Fig.|9]for an 
illustration.) The behavior of these roots is as follows. For > 2fW^' cubic has no non-negative real roots. 
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Figure 9: Cubic f{u) = 2Mv? — v? + in denominator of Eq. ((T3]l . The red curve has (? = ffM^' while the 
green has = 0. 

At (? = TyfKJ^ ^ double root at n = Jfrr (and a third at ti = which is irrelevant because we consider 



only non-negative values of = l/r). This double root, which causes the integrand to diverge, corresponds to 
the pure angular orbit at n = 3^. As we further reduce c^, the double root aiu = separates into two, the 
larger increasing toward u = (the horizon), and the smaller decreasing toward u = (asymptotic infinity), 
which they reach when (? descends to its smallest possible value 0. So the task of our algorithm is to shoot null 
geodesies from E\ in different directions (different values of c^), trying to hit 7 ((/72)/ all the time being careful 
to keep (? large enough that the roots of the cubic do not fall between either u\ or U2, as doing so would cause 
the integrand to become infinite or imaginary. 

Our algorithm begins by setting = if either event lines behind the horizon, or (? = i^fM^ otherwise. It 
will later adjust the value of (? , in its attempt to locate a null geodesic which travels from E\ to 7. For each 
value of (? , it first checks to see if either u\ or U2 lies between the non-negative real roots (if any). If so, it adjusts 
(? upward (by using a linear extrapolation of j{u) at the roots), until neither u\ or U2 lies between them. 

Now that we have an integrand which is real and finite in the entire domain [ui, ii2]/ the code integrates 
Eq. ((T3|) numerically using the composite Simpson's rule, with n = 512 subintervals. If either the numerically 
evaluated integral is within A = 5e of ip2, where e is double precision machine epsilon (srf 2.22045 x 10^^^ on 
the machine on which the results of section |3] were generated), or the two most recently chosen values of 
are within 5e of each other (see below), we then (tentatively) decide that this null geodesic arrives at 7. If this 
condition fails, we choose a second value for and repeat. The second value is chosen to be the the first + .03 
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if the integral overshoots (is greater than) ip2, or the first - .005 if it undershoots. If the second guess of (? also 
misses Lp2, then subsequent values are chosen by a linear interpolation/extrapolation from the two previous 
guesses. This algorithm (with the additional features described below) converges for all pairs of events we 
have encountered in our simulations. 

There is a special situation which can arise (as in some of the examples of Subsection 13. l)) , in which there 
are no causal curves from E\ to 7. This occurs when E2 lies behind the horizon, and the integral with = 
undershoots Lp2- This means that every future directed causal curve from E\ falls into the singularity before 
reaching 7, so the events must be unrelated. 

When (? is small enough that the domain of integration touches a root, the integral diverges. Often the 
'target value' of 992 requires a (? which is is very close to this singular value. We find that a convenient way to 
handle this situation numerically is to detect when we manage to find a valid value of (? , which is large enough 
for u\ and U2 to escape the roots, and yet small enough to yield an integral which exceeds Lp2- Once we find 
this value of (? = c^jj^, then we know that the value of (? we seek is greater than this. Thus, in the course of 
the above iteration, which uses linear interpolation/ extrapolation to select subsequent values of (? , if a value 
is selected which is smaller than (? = c^^^j^, then we instead choose the mean of = c^j^ and the previous c^. 
(Furthermore in subsequent iterations, if we find a yet larger value of for which the integral exceeds ^p2, then 
we use this as the new c^j^.) 

Once the above loop converges, so that we have a value of for which the integral of Eq. ([131 ) yields 1^2, 
we then check that the numerical approximation to the integral is sufficiently accurate. The check is simple: 
we compute the numerical approximation to the integral again at four times the resolution (/92(4n) (4n subin- 
tervals), and subtract that value from the 'P2{'n) using n subintervals. If the difference is greater than Srj, where 
•q is the larger of tp2 (4n) — (/?2 and 5e, then we double n and repeat the above iteration. (Though we stop the 
iteration if the difference |(/?2(4n) — <f2{'n)\ ever increases from that for the previous n.) 

Now that we have an accurate value of c^, which yields a null geodesic which hits 7, we integrate Eq. ((15)) 
to get the elapsed time along the geodesic, and thus can determine if the two events are related by comparing 
this elapsed time with t2 —ti. 
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